Anomalous synchronization threshold in coupled logistic maps 
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We consider regular lattices of coupled chaotic maps. Depending on lattice size, there may exist a 
window in parameter space where complete synchronization is eventually attained after a transient 
regime. Close outside this window, an intermittent transition to synchronization occurs. While 
asymptotic transversal Lyapunov exponents allow to determine the synchronization threshold, the 
distribution of finite-time Lyapunov exponents, in the vicinity of the critical frontier, is expected 
to provide relevant information on phenomena such as intermittency. In this work we scrutinize 
the distribution of finite-time exponents when the local dynamics is ruled by the logistic map 
1 h> 4i(l - i). We obtain a theoretical estimate for the distribution of finite-time exponents, that 
is markedly non-Gaussian. The existence of correlations, that spoil the central limit approximation, 
is shown to modify the typical intermittent bursting behavior. The present scenario could apply to 
a wider class of systems with different local dynamics and coupling schemes. 



PACS numbers: 05.45.Ra,05.45.-a,05.45.Xt 

I. INTRODUCTION 

Coupled map lattices (CMLs), dynamical systems with 
discrete space and time, are being intensively investi- 
gated since the early 80's as models of many spatiotempo- 
ral phenomena occurring in a wide variety of systems [l| . 
One of such phenomena is synchronization and, in partic- 
ular, amongst the various kinds of synchronized behavior, 
complete synchronization (CS) 2] occurring in regular ar- 
rays of coupled chaotic systems. Intermittency plays an 
important role in the destabilization of the synchronized 
stated . Intermittent transitions between laminar, quies- 
cent states, and irregular bursting have been investigated 
in many dynamical systems, their scaling properties be- 
ing determined 0, IS @- Theoretical derivations of these 
scalings usually assume Gaussian fluctuations for describ- 
ing the irregular bursts between consecutive laminar re- 
gions. Although this may be a suitable approximation 
for many systems, there are cases where deviations from 
Gaussianity can be observed, specially in the vicinity of 
the synchronization threshold. It is our purpose to ex- 
hibit such anomalous behavior. 

Let us consider, as paradigmatic example of a spatially 
extended system, periodic chains of N one-dimensional 
chaotic maps x i— ► f(x) evolving through a distance- 
dependent diffusive coupling: 

N' 

x\% = (l- £ )/(^)+£^S( r ) (f(xt r) ) + f(4 i+T) )) , 

^ r=l 

(1) 

(i) 

where x t represents the state variable for the site i 
(i = 1, 2, N) at time t, e > is the coupling strength, 



B(r) an arbitrary function of r and n — 1 B(r) a 

normalization factor, with N' = (N — l)/2 for odd N. 

CS takes place when the dynamical variables that de- 
fine the state of each map adopt the same value for all 
the coupled maps at each time step t, i.e., = x[ 2 ' = 
... = Xt = x[*\ It can be easily verified that this 
state is solution of Eq. (Q. The question is whether 
it is stable or not with respect to small perturbations 
in directions transversal to the CS state in CML phase 
space. A criterion for stability can be drawn from the 
asymptotic Lyapunov spectrum (LS). As far as the CS 
state lies along the direction associated to the largest ex- 
ponent, the CS state will be transversely stable if the 
(N — I) remaining exponents are negative. Therefore, if 
there is a single attractor, negativity of the second largest 
Lyapunov exponent implies that the CS state is asymp- 
totically attained^. Depending thus on system size and 
on the parameters defining both the chaotic map and 
B(r), there may exist an interval of values of the coupling 
strength e for which the CS state is asymptotically sta- 
ble. Explicit results have been shown before, for instance, 
for algebraically decaying interactions, i.e., B(r) = l/r a , 
with < a|3,|9( and for uniform interactions with a cut- 
off distance /?, with 1 < (3 < jV'pjJ. In those cases, for a 
given system size, one finds a synchronization domain in 
parameter space. 

Inside this domain, CS eventually occurs after a tran- 
sient whose typical duration diverges as one approaches 
the critical frontier. Outside the synchronization domain, 
there may occur intermittent behavior. A characteriza- 
tion of phenomena associated to a blowout bifurcation, 
such as intermittency [Tl| . can be performed through the 
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analysis of the distribution of largest transversal finite- 
time Lyapunov exponents (LTFEs). We will show that 
for Ulam maps coupled through schemes of the form , 
the probability density function (PDF) of LTFEs devi- 
ates from a Gaussian law. This will be shown to have 
important consequences at criticality. Let us remark that 
analytical results are specially relevant in this field due to 
the very nature of the phenomena involved, sensitive to 
the unavoidable finite precision of numerical simulations 




II. FINITE-TIME LYAPUNOV SPECTRUM 

The dynamics of tangent vectors £ = (faW, Sx^ 2 \ . . . , 
5x( N ^) T is obtained by differentiation of the evolution 
equations Q). In order to obtain finite-time exponents, 
we proceed analogously to what we have done before 
for the calculation of asymptotic exponents -8, 9J. For 
self-containedness, let us review the steps. The tangent 
dynamics is given by = Tj£j, where the Jacobian 
matrix Tj is 

T t = ((l-eJ + ^BjDi, (2) 

with the matrices Dj and B defined, respectively, by 
Dj k = f{x\ j) )6 jk and B jk = B{r jk ){l - 6 jk ) , being 
Tj k = miiii e z\j — k+lN\. Once the initial conditions have 
been specified, the LS of finite-time Lyapunov exponents 
(FTEs) {Afc(n)}, calculated over a time interval of length 
n, is extracted from the evolution of the initial tangent 
vector £ : £n = T„£ , where T n = T n _i . . . TiT . The 
FTEs are obtained as X k (n) = In An , for fc = 1 , . . . , TV, 
where {A„ } are the eigenvalues of A n = (TjT„ ) ^ |l3| . 

In CS states, the dynamical variables of all maps coin- 
cide at each time step i, i.e, x^ = x^ = . . . = = 
x ( * } . In this case, D, = /'(if'jUv, thus, Tj = /'(^* } )B 
and TT% = (S$l£[f'(xW)] 3 )B 2i . Therefore, one ar- 
rives at the following expression for the spectrum of N 
finite-time Lyapunov exponents, over a time interval of 
length n, in CS states: 

-. n — 1 

Afc(n) - -y2in\f'(xV)\+c k , forfc=l,...,JV, (3) 

i=o 

with Cfc = In 1 1 — e + eb k /rj\, where b k , the eigenval- 
ues of the interaction matrix B, are given by b k — 
2J2^ =1 B(m)cos(2Trkm/N), for odd N. Here x { * ] = 

P{xq ) is the ith iterate of the initial condition Xq , the 
same for the N maps, since we are dealing with CS states. 

In the asymptotic case n — * oo, and assuming ergodic- 
ity, the first term in the left-hand side of Eq. J2J), which 
represents a time-average, can be substituted by an aver- 
age over the single-map attractor. In such case one gets 
the asymptotic LS 

Afc = Xu + Cfc, (4) 



where Xu — (In is the Lyapunov exponent of 

the uncoupled map. Notice that the parameters that 
define the particular uncoupled map affect only Xu, while 
{cfc} are determined by the particular cyclic dependence 
on distance in the regular coupling scheme. 

As discussed in the Introduction, the asymptotic LS 
provides a criterion for synchronization, namely the neg- 
ativity of second largest (or largest transversal) asymp- 
totic exponent Aj_- For the particular interaction B(r) = 
1 /r a , that we will consider in numerical simulations, the 
above condition leads to|l,|3 

1 - c~ Xu 1 + c~ Xu 

III. DISTRIBUTION OF FINITE-TIME 
LYAPUNOV EXPONENTS 

Unlike infinite-time exponents, local exponents 
strongly depend on the initial conditions. Starting from 
random initial conditions, the fluctuations in the values 
of Afc(n) arise from the summation in @ only, since 
{cfc} are constant. As a consequence, all the X k (n) of 
the spectrum will have a PDF with the same shape, 
differing only in the mean value (X k (n)) — Xu+c k , which 
coincides with the asymptotic exponent A^. In particular 
the variance of finite-time Lyapunov exponents (FTEs) 
in the CS states does not depend on c k (since it is an 
additive constant). Therefore, it is the same for all the 
exponents of the spectrum, as expected also from the 
fact that the source of the fluctuations is unique. Along 
a trajectory, fluctuations shift the finite-time Lyapunov 
spectrum as a whole. A further consequence is that the 
variance of the local exponents in CS states does not 
depend on the lattice parameters embodied in {c k }, but 
only on the local features of the individual map. In fact, 
in CS states, all maps evolve with the dynamics of an 
uncoupled map. Being so, the PDFs of FTEs can be 
straightforwardly obtained from the PDF of the local 
map. Hence, let us omit by the moment the index k. 
Our discussion will be valid for any X k (n), in particular 
for X±(n). 

For the Ulam map, x i— > 4x(l — x), a smooth approx- 
imate expression for the PDF of FTEs is available [l4| . 
namely 

In 

P(X(n)) ~ — ln(coth |n[A(n) - A]/2|) , (6) 

valid for large enough n. For X(n) < A, expression © 
is exact, however, for A(n) > A, the exact distribution 
presents a complex structure with 2™" 1 spikes that get 
narrower and accumulate close to the mean value with 
increasing n0, Therefore, in the latter interval, 

expression © constitutes a smooth approximation, such 
that the sharp spikes have been trimmed by finite size 
bins. Moreover, the exact PDF is non-null only for X(n) — 
X < In 2. 
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Even the smooth distribution © is markedly different 
from Gaussian. It is divergent at A(n) = A and falls off 
with exponential tails. It is noteworthy that the variance 
decays anomalously as 1 /n 2 0, 0] 

instead of the usual 1/n decay. Also notice that the PDFs 
© for different values of n would collapse into a single 
shape via rescaling by n. 
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FIG. 1: Comparison between numerical and theoretical PDFs. 
They correspond to the largest transversal time-20 exponent 
of N — 21 coupled Ulam maps, with algebraically decay- 
ing interactions, for e = 0.8 and different values of a (above 
a c — 0.867 stability of the CS state is lost). Numerical PDFs 
are represented by symbols. Solid lines associated to full sym- 
bols correspond to the theoretical prediction given by Eq. 0, 
while those associated to hollow symbols correspond to Gaus- 
sian fittings. Inset: semi-log representation to exhibit the 
tails. 

Fig^exhibits numerical PDFs for CMLs together with 
the analytical prediction given by Eq. ©. Numerical 
PDFs were built by choosing 10 4 initial conditions and 
computing the second eigenvalue of the matrix A„ after 
a transient. Eq. JSJ), obtained for uncoupled maps, is in 
excellent agreement with numerical results for LTFEs in 
the coherent states of CMLs, as expected. It is worth 
noting that, although expression © was derived for CS 
states, it remains still a good approximation even close 
outside the synchronization domain. This means that the 
correlations that lead to violation of the central limit the- 
orem persist in that region. However, far enough from the 
threshold, the terms that contribute to FTEs become un- 
correlated by the "bath" of coupled chaotic maps, and the 
central limit theorem holds, leading to Gaussian shapes 
(see Fig.Q}. Notice that the PDF for a = 1.0, although 
Gaussian in the central part, stills falls at right with a 
exponential tail. A detailed analysis of correlations giv- 
ing rise to JfJJl can be found in Refs. 0,0]- Basically, 



the autocorrelation function C(n) of one-step FTEs is 
C(n) = — 7r 2 2~™/24, that is, it decays exponentially with 
time n (with a characteristic time equal to the inverse 
of Xu = In 2). However, since successive one-step ex- 
ponents are dependent through the deterministic logistic 
mapping, non-Gaussian PDFs of FTEs arise. 



IV. CONSEQUENCES OF FLUCTUATING 
EXPONENTS 

A. Subcritical regime 

For parameter values belonging to the synchronization 
domain, that is, for subcritical a, being all other param- 
eters fixed, the system eventually converges to the CS 
state. In fact, asymptotically, CS states are stable since 
the distribution of LTFEs collapses to a Dirac delta func- 
tion centered at Aj_, which is negative in that domain. 
The relaxation to a CS state can be measured, for in- 
stance, by means of either the distance to the SM, defined 

through d(t) = — (x t )) 2 , or the order param- 

eter R(t) = |X; i exp(27rxf ) )|/A. Since, for small devia- 
tions from the SM, both quantities are related through 
d 2 ~ (1 — i? 2 ) /2, we will exhibit the time evolution of 
d 2 only. After a very brief transient, the decay to the 
CS state is exponential with a characteristic time given 
by t c = 1/| A j. | (see Fig. that diverges at the criti- 
cal frontier. For the power-law interaction, A^ scales as 
|Aj_| ~ \ot — Oi c \ and |Ax| ~ \e — e c |, at the critical point. 
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FIG. 2: Relaxation to the synchronization manifold. Time 
series of d 2 for the same CML of Fig. Qwith a — 0.8 (sub- 
critical). The dashed line corresponds to the exponential law 
indicated in the figure, for comparison. 

The fact that the distribution of LFTEs spreads over 
negative and positive values (Fig. QJ, implies that the 
exponents, computed over finite-time segments of a tra- 
jectory, fluctuate around zero. On one hand, as one 
approaches the frontier subcritically, the mean value of 
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the distribution shifts to zero from negative values. On 
the other hand, as one follows a trajectory for a longer 
time interval, the PDF of LTFEs concentrates around 
the mean. Then, there may be segments of trajectory 
in which the lattice is repelled from the synchronization 
manifold (SM). But, in average, it is attracted exponen- 
tially fast. Due to the finite precision of computer cal- 
culations, the distance to the SM saturates (see Fig- EJ - 
Intrinsic noise, due to numerical truncation, may drive 
the state of the system slightly away from the saturation 
level. However, each time this happens, the distance de- 
cays, again exponentially fast, to its lower bound. 

Numerical results in the subcritical regime, illus- 
trated in Fig. El can be suitably described by d(t) = 
d exp[Aj_(£) t]. Then, from Eq. ®, we have 



d(t) = d exp[\±t + ^2c(n)] 



(8) 



where we have split the fluctuating component £ = 
A_l(1) — A_l, corresponding to successive (time-correlated) 
centered one-step LTFEs. The evolution of the distance 
to the SM can also be modeled by a ltd multiplicative 
stochastic differential equation for d, 



d = \ ± d + C(t)d, 



(9) 



where £(t) is a colored non-Gaussian noise and time is 
continuous. 

As shown in Sect. IIIII in particular, the distribution 
of Aj_(l) in synchronized states follows that of one-step 
FTEs in the uncoupled map. An exact expression can be 
straightforwardly derived from the invariant measure of 
the Ulam attractor^5|. Then, one has 



no 



i 



7T V4c- 2< ; - 1 



(10) 



with ( S (—oo, In 2). Notice that £, with zero-mean, is 
bounded from above. This explains the upper bound of 
the fluctuations superposed to the exponential decay in 
Fig. J5J. In fact, the fluctuations around the exponential 
envelope can be fully described by the statistics of time- 
one exponents. The time series and distribution of (, are 
presented in Fig. [3| for supercritical a, showing that even 
close outside the synchronization domain, the distribu- 
tion of £ follows that of the uncoupled map. 

Breakdown of shadowability might occur 0, as ex- 
ponents fluctuating about zero are a signature of un- 
stable dimension variability. Taking into account that 
the PDF is non-null for Aj_(n) < A i + ln2[Tl H^. 
then if Aj_(a,e, N) < —In 2, the finite-time exponents 
are negative for almost any initial condition. Thus, 
this point would correspond to the onset of shadowa- 
bility. Whereas, for — In 2 < X±, although the mean 
of the distribution may be negative, there is always a 
non-null fraction / of positive exponents given by / = 
J °° dX±(n)P(X±(n)), pointing to the possibility of loss 
of shadowability of numerical trajectories^). Because / 
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FIG. 3: Time evolution of £ = Ax(l) — Aj_ (A) for the same 
CML of Fig. with a = 0.9 (supercritical). Normalized his- 
togram of £ (B). Inset: semi-log representation to exhibit the 
exponential tail. Full lines correspond to Eq. HUB . 



grows from zero with a very small slope, since the posi- 
tive tail of the PDF is approximately exponential, then 
the onset may appear shifted towards the threshold in 
numerical computations fl9| . 



B. Supercritical regime 

For supercritical a, that is outside the synchronization 
domain, CS states are not asymptotically stable, because 
Aj_ > 0. Close to the threshold (up to a ~ 1 for e = 0.8), 
the numerical estimate for the second largest transversal 
exponent fairly coincides with Aj_ (calculated for synchro- 
nized states). Moreover, close enough to the boundary, 
correlated bursts away from the SM occur (see Fig. 0}. 
Although this figure exhibits a time series up to t = 2000, 
the same features are observed for longer runs (performed 
typically up to t = 10 7 ). 

The intermittent behavior can also be understood 
in terms of the distribution of the finite-time largest 
transversal exponent Aj_(n), that close to the thresh- 
old can be described by Eq. |JBJ|. Although the average 
(Aj_(n)) = A^ is positive, there is a non-null probability 
that the LTFE be negative, thus leading to intermittent 
behavior at finite times. The distribution of LTFEs indi- 
cates that there are time intervals during which the tra- 
jectories are either attracted to or repelled from the SM. 
Despite the average duration of the synchronized time in- 
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tervals increases when approaching the critical frontier, 
synchronization is not attained as a final stable state. 
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FIG. 4: Time series of d 2 for the same CML of Fig. □ with 
a ~ 0.9 (supercritical). Inset: semi-log representation of the 
same data. The dashed line corresponds to the exponential 
law indicated in the figure, for comparison. 



where £ is a fluctuating quantity with zero-mean and 
variance af. If fluctuations were white Gaussian ones, 
then, from the stationary solution of its associated 
Fokker-Planck equation, the distribution of y should in- 
crease following the law P(y) ~ exp(hy), with h de- 
fined as above. Multiplicative corrections to Eq. Qlip. 
to model the upper-bounded behavior of d, will affect 
the shape of P(y) mainly above the average value, and 
are not expected to affect the small d behavior. The main 
point is that the stochastic process {£} is not white nor 
Gaussian, as can be neatly observed in Fig. [21 what may 
explain the deviation from the exp(/iy) law. 




Intermittent bursts grow exponentially with a charac- 
teristic time given by 1/X± (see inset of Fig. 0J). There- 
fore, the characteristic time has the same scaling laws as 
the subcritical characteristic time discussed above. The 
distance to the SM fluctuates around a reference level 
((d) ~ 0.037 for the parameters in Fig. 0} that increases 
with a. Nonlinearities in Eq. © cause saturation of the 
growth experienced by the distance to the SM due to 
the positiveness of A_l, and keep the distances within a 
bounded interval. Close to the threshold, the average 
distance increases linearly with Aj_ ~ \a — a c \. Also, for 
increasing a, the correlated bursts become more frequent 
(hence, its duration becomes shorter) such that far from 
the frontier fluctuations become uncorrelated and the in- 
termittent clustering effect disappears. Moreover, in that 
region, the numerical estimate of the largest transversal 
exponent significantly deviates from the one calculated 
for synchronized states and new features, out of the scope 
of the present work, occur. Therefore, we will focus on 
the near vicinity of the threshold. 

The histogram of logarithmic distances y = In d is pre- 
sented in Fig. [S] The distribution initially grows approx- 
imately as exp(az/) and above the maximum value falls 
off faster than exponentially. The coefficient of the ex- 
ponential argument is a ~ 1, at variance, with the value 
a = 2X±/a 2 = h (hyperbolicity exponent) 6], where a 2 is 
the variance of time-1 exponents, derived under assump- 
tion of Gaussian fluctuations. In fact, in the continuum 
approximation, the evolution of y = lnd follows, at first 
order, the Langevin equation 

y = \ ± + C(t) (11) 



FIG. 5: Distribution of the logarithmic distance y — Ind, 
for different values of a. Inset: exponent a, resulting from 
the best exponential fit to the initial increasing regime, and 
hyperbolicity exponent h and as a function of a. 

A universal result is that, at the onset of intermit- 
tency, the distribution of laminar phases (inter-burst in- 
tervals) decays as a power-law0, meaning the presence 
of lengths of arbitrarily large size. In particular the size 
of the average plateau diverges. Moving far from the 
onset, the tail of the distribution of laminar phases is 
gradually dominated by an exponential decay. In gen- 
eral, the distribution of lamellar phases can be obtained 
by solving a first-return problem. For the power-law de- 
cay, the exponent j3 = 3/2 was found to be universal, as 
far as the central limit approximation holds 4]. However, 
as we have seen, in our case, very close to the threshold, 
correlations persist and the distribution of finite-time ex- 
ponents, even, and specially, in the central part (where 
it is divergent), deviates from the Gaussian approxima- 
tion. Therefore, deviations from the 3/2 power law are 
expected. We measured inter-burst sizes, that is, the 
length of time segments during which the distance d re- 
mains below a threshold value d . Numerical distribu- 
tions of inter-burst sizes, for values of parameter a close 
to the synchronization threshold, are displayed in Fig. [S] 
In general, a rapid decay is observed for very small r 
(not exhibited). Since a power-law behavior is expected 
for sufficiently large r, histogram heights for small val- 
ues of t are not exhibited in the figure, for the sake of 
clearness. One sees that the PDF of inter-burst intervals 
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follows a power-law with j3 ~ 1, at neat variance with the 
value (3 = 3/2. A cross-over to a asymptotic exponen- 
tial regime is always observed. We have already seen in 
Fig- CI that, below a ~ 1, correlations persist. This may 
explain the power-law exponent observed in Fig. [5] 
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V. SUMMARY 

We have presented analytical results for the PDF of 
LTFEs in CS states of coupled Ulam maps. They were 
confirmed by numerical experiments performed for CMLs 
with interactions decaying with distance as a power law. 
The knowledge of the statistical properties of finite-time 
Lyapunov exponents allows to understand the anoma- 
lies in bursting behavior close to the SM. As a conse- 
quence, universal laws derived under assumption of Gaus- 
sian fluctuations do not hold generically. 

Our results were obtained for local Ulam maps, how- 
ever, the observed features may occur in a wider class of 
extended systems, as long as other deterministic chaotic 
or stochastic processes present similar deviations from 
Gaussianitv[l5|. 



FIG. 6: Distribution of inter-burst times for the same CML 
of Fig. with different (supercritical) values of a indicated 
on the graph. In all cases, d — 2(d), but we verified that 
decay laws do not substantially depend on the choice of the 
threshold d . Full lines, corresponding to the indicated power- 
laws, were drawn for comparison. 



Acknowledgments 

This work was partially supported by Brazilian agen- 
cies CNPq and Fundacao Araucaria. 



[1] K. Kaneko, in Theory and Applications of Coupled Map 
Lattices, edited by K. Kaneko (Wiley, Chichester, 1993). 

[2] A. Pikovsky, M. Rosenblum, J. Kurths and S. Strogatz, 
Synchronization: A universal concept in nonlinear sci- 
ences, (CUP, Cambridge 2001); S. Boccaletti, J. Kurths, 
G. Osipov, D.L. Valladares and, C.S. Zhou, Phys. Rep. 
366, 1 (2002). 

[3] N. Piatt, E.A. Spiegel and C. Tresser, Phys. Rev. Lett. 
70, 279 (1993) 

[4] J.F. Heagy, N. Piatt and S.M. Hammel, Phys. Rev. E 49, 
1140 (1994); H.L. Yang and E.J. Ding, Phys. Rev. E 54, 
1361 (1996). 

[5] A. Cenys, A. Namajunas, A. Tamasevicius and T. Sch- 

eneider, Phys. Lett. A 213, 259 (1996); M. Sauer and F. 

Kaiser, Phys. Lett. A 243, 38 (1998). 
[6] S.C. Venkataramani, T.M. Antonsen Jr., E. Ott and J.C. 

Sommerer, Physica D 96, 66 (1996). 
[7] P.M. Gade and C.-K. Hu, Phys. Rev. E 60, 4966 (1999); 

P.M. Gade, HA. Cerdeira and R. Ramaswamy, Phys. 

Rev. E 52, 2478 (1995). 
[8] C. Anteneodo, S.E. de S. Pinto, A.M. Batista and R.L. 

Viana, Phys. Rev. E 68, 045202(R) (2003); ibid. 69, 

029904(E) (2003). 



[9] C. Anteneodo, A.M. Batista and R.L. Viana, Phys. Lett. 

A 326, 227 (2004). 
[10] P.G. Lind, J. Corte-Real, J.A.C. Gallas, Phys. Rev. E 

69, 026209 (2004). 
[11] M. Ding and W. Yang, Phys. Rev. E 56, 4009 (1997). 
[12] C. Grebogi, S.M. Hammel, JA. Yorke and T. Sauer, 

Phys. Rev. Lett. 65, 1527 (1990); C. Zhou and J. Kurths, 

Phys. Rev. Lett. 88, 230602 (2002). 
[13] J.-P Eckmann and D. Ruelle, Rev. Mod. Phys. 57, 617 

(1985). 

[14] C. Anteneodo, Phys. Rev. E 69, 016207 (2004). 
[15] A. Prasad and R. Ramaswamy, Phys. Rev. E 60, 2761 
(1999). 

[16] J. Theiler and LA. Smith, Phys. Rev. E 51, 3738 (1995). 
[17] S. Dawson, C. Grebogi, T. Sauer, and J. A. Yorke, Phys. 

Rev. Lett. 73, 1927 (1994); TD. Sauer, Phys. Rev. E 65, 

036220 (2002). 

[18] R.L. Viana, S.E. de S. Pinto, J.R.R. Barbosa and C. 

Grebogi, Int. J. Bifurcat. Chaos 13, 3235 (2003). 
[19] R. L. Viana, C. Grebogi, S. E. de Souza Pinto, S. R. 

Lopes, A. M. Batista, and J. Kurths, Phys. Rev. E 68 

(2003). 



